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We consider the one-dimensional motion of a particle randomly accelerated by 
Gaussian white noise on the line segment < x < 1. The reflections of the particle 
from the boundaries at x = and 1 are inelastic, with velocities just after and 
before reflection related hy Vf = —rvi. Cornell et al. have predicted that the 
particle undergoes inelastic collapse for r < Tc = e^'^l^ = 0.163, coming to rest at 
the boundary after an infinite number of collisions in a finite time and remaining 
there. This has been questioned by Florencio et al. and Anton on the basis of 
simulations. We have solved the Fokker-Planck equation satisfied by the equilibrium 
distribution function P(x, v) with a combination of exact analytical and numerical 
methods. Throughout the interval < r < 1, P{x,v) remains extended, as opposed 
to collapsed. There is no transition in which P{x, v) collapses onto the boundaries. 
However, for r < r^ the equilibrium boundary collision rate is infinite, as predicted 



by Cornell et al., and all moments g > of the velocity just after refiection 
from the boundary vanish. 
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I. INTRODUCTION 



Consider a particle randomly accelerated on the line segment < x < 1 according to 

(fx 



^ = v{t). {v{tHt')) = 26{t-t'), (1) 



where ri{t) is uncorrelated white noise with zero mean. If the collisions of the particle with 
the boundaries at x = and 1 are elastic, the mean square velocity increases according to 

{v{tf) = {vm+2t, (2) 

just as in the absence of boundaries. 

In this paper we assume that the boundary collisions of the randomly accelerated particle 
are inelastic. The velocities just after and before reflection satisfy 

Vf = -rvi, (3) 

where r is the coefficient of restitution. This simple model is of interest in connection with 
the statistics of driven granular media, where particles tend to cluster, due to inelastic 
collisions, even though no attractive forces are present. The model was studied by Cornell, 
Swift, and Bray (CSB) who argued that the particle undergoes "inelastic collapse," i.e. 
makes an infinite number of collisions in a finite time, comes to rest at the boundary, and 
remains there, if the coefficient of restitution r is less than the critical value 

r^ = e-'^/^ = 0.163... (4) 

The prediction of inelastic collapse was questioned by Florencio et al. P, who carried 
out simulations and found that the particle did not adhere to the boundary for any r. Anton 
reported that his simulations are consistent with an infinite collision rate for r < Tc but 
also incompatible with localization of the particle at the boundary. 

According to Eqs. (j21) and the kinetic energy of the randomly accelerated particle 
increases in between boundary collisions but decreases, for r < 1, in the collisions. Eventu- 
ally an equilibrium is reached. Burkhardt, Franklin and Gawronski (BFG) analyzed the 
equilibrium distribution P{x,v) for the position and velocity of the particle for rc < r < 1. 
This function satisfies the steady-state Fokker-Planck equation 

d 92 \ 
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with the boundary conditions 

P{x,v) = P{l-x,-v), (6) 
P{0,-v) = r^P{0,rv) , v>0, (7) 

corresponding to reflection symmetry and conservation of probability, respectively. In par- 
ticular, the second boundary condition insures that the incident and reflected probability 
currents at the boundary have equal magnitude 

poo rco 

1=1 dvvP{0,-v)= / dvvP{0,v). (8) 
Jo Jo 

Making use of an exact Green's function solution of Eqs. ©-©i BFG found that the 
boundary collision rate /, defined by Eq. (jHJ, diverges as r approaches from above and 
that P{x,v) is extended, as opposed to collapsed, at r = Vc- In this approach P{x,v) is 
obtained as the difference of two integrals, both of which diverge for r < re- This was 
noted by BFG, who, however, incorrectly concluded that the solution to the Fokker-Planck 
equation breaks down for r < r^. 

In this paper the calculation of P{x, v) is extended to r < r^. In Section II the approach 
of BFG is reviewed. The divergences, for r < r^, of the two integrals which determine P{x, v) 
are shown to cancel, leaving a finite result. Throughout the entire interval < r < 1, P{x, v) 
varies smoothly and analytically with r. There is no transition in which P{x, v) collapses 
onto the boundaries. 

In Section III the equilibrium boundary collision rate is calculated from the results of 
Section II. The collision rate is finite for r > and infinite for r < r^, as predicted by CSB. 
All the equilibrium moments |f |'^, q > of the velocity just after striking the boundary ^ 
vanish for r < Vc- 

Our conclusions are summarized in Section IV, and some earlier results on inelastic col- 
lapse are reexamined. 



II. SOLUTION OF THE FOKKER-PLANCK EQUATION 

We begin with a brief review of the approach P| of BFG. Generalizing earlier work of 



Masoliver and Porra 6\, they showed that the Fokker-Planck equation ©, with reflection 
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symmetry ©, has the exact solution 

P{x, v) = 

for f > in terms of the Green's function 

yl/2^1/2 



duuG{x, V, u)P{0, u) 



(9) 



G{x, V, u) 



{»;3+u3)/9x 



3x 



/2y3/2^3/2X 



1 



3i/3r(2 
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dy- 



/3 



V 9x 



where 



R{y,u) 



^l/2g-«V% 



35/6r(i)r(|) yV6(i_^)i/6 



T76 1^1 



1 5 u^{l-y) 



6'6' 



91/ 



(10) 



(11) 



and iFi{a,b,c) is a standard confluent hypergeometric function 

To calculate f ) from Eq. 0, one must first determine the unknown function P(0, u) 
on the right hand side. Setting a; = 1 in Eqs. ©-(HH) and using r^P{0,rv) = P{l,v), as 
follows from Eqs. (0) and ((Tj), leads to the integral equation 



POO 

r^P{0,rv)= duuG{l,v,u)P{0,u) 
Jo 



for P(0, v), where 

V, u) 



.^1/2^1/2 g-(„3+«3)/9 
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+ 61F2 1; 



5 7 f^M^ 



'6' 6' 81 



(12) 



(13) 



and iF2{a; b, c; z) is a generalized hypergeometric function The quantity ^^(l, v, u) is of 
interest in its own right. As discussed in the Appendix, it generalizes McKean's result 
for the velocity distribution at first return to the boundary from the half line a; > to the 
line segment < x < 1. 

BFG showed [4] that the asymptotic form of P(0, v) for small and large v is determined 
by the first and second terms, respectively, of the kernel G{l,v,u) in Eqs. (jl2j) and (jl3j) and 
given by 

, ^ , 



P{0,v) 



(14) 



where 



2 sin 



2/3+1 



■ TC 



l/(/3-2) 



9^3 
1 — r'^ 



(15) 
(16) 
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Note the non-Maxwellian velocity distribution. As r decreases, the boundary colhsions 
become more inelastic, and the probability of finding the particle near the boundary with a 
small velocity increases. This is seen in the monotonic increase of the exponent /5(r) from 
to |, as r decreases from 1 to 0. The characteristic velocity fch('") also decreases with 
decreasing r. 

The asymptotic forms ()14j ) - (fT^ are smooth analytic functions of r throughout the interval 
< r < 1. There is no singular behavior at Vc- In particular, on expanding the right side 
of Eq. (fT3j) about P = 2, one sees that /9(r) is a nonsingular function of r, with /?(rc) = 2. 
To connect the asymptotic forms ()14|l - (fT?)|l of P(0,f)for small and large v, we have solved 
integral equation (jl2|) by numerical iteration, as in ^. As noted by BFG [J], the integral 
equation appears to have a well-defined solution for < r < 1, i.e. < /? < 5/2, with no 
special behavior at r^. Numerical results for several values of r above and below Tc = 0.163, 
are shown in Fig. 1. The slopes of the curves, for small v, depend on r in accordance 
with the asymptotic form (|T4j) . |T5|) . There is no qualitative difference above and below Tc- 
Presumably, P{0,v), like its exact asymptotic forms (fT^ - |TH| for small and large v, is an 
analytic function of r throughout the interval < r < 1. 

Once P(0,f) has been determined, P{x,v) may be obtained by integration. According 
to Eqs. (P)-(fTT|). P{x,v) is the sum of two integrals over P{0,u), corresponding to the two 
terms in the Green's function G{x,v,u) in Eq. (fTUI) . Both integrals diverge at the lower 
limit for r < Tc, as follows from the asymptotic^ form (|T4j) . ()15|1 of P{0,u) for small u, with 
/3(r) > 2 for r < re- This was noted by BFG [4], who, however, incorrectly concluded that 
the solution to the Fokker-Planck equation breaks down for r < r^. The divergences cancel, 
leaving a finite result, as may be seen by integrating with a low u cutoff and sending the cutoff 
to zero after adding the two integrals. No cutoff is needed if the two integrands are added 
before integrating over u. From Eq. (fTUI) it is straightforward to show that G{x, v, u) ~ u^f^ 
in the small u limit Thus the integral in Eq. (jHl) behaves as j^duu^^"^'^ for small u. 
Since < (3{r) < | for < r < 1, there are no convergence problems at the lower limit of 
the integral. Throughout the interval < r < 1, P{x, v) is a smooth well-defined function of 
r, presumably analytic in r, and does not collapse onto the boundaries at x = and a; = 1. 

We have also considered the probability density 

/oo poo 
dvP{x,v)= dv [P{x,v) + P{1- x,v)] (17) 
oo Jo 
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for the position of the particle. From Eqs. Q-dJ and ()17p. one finds that the leading 
singular contribution to P{x) for a; ^ is determined by the asymptotic form P(0, v) ~ 
Av~^ for V ^ in Eq. ^T^i and given by 

P^i^^ix) ^ Bx^'-"^/' , x^O, (18) 
2vr r(^) 

3(4/3+5)/6 ^) r(f )r(^) ■ ^ ^ 

For < /3 < 1, i.e. | < r < 1, the leading singular contribution to P{x) = P(l — x) in 
Eq. ()18|) vanishes as x approaches or 1, and -P(O) is finite and nonzero. For /5 > 1 or 
r < i, P{x) diverges according to Eq. (fTHj) as x approaches or 1. Since the divergence is 
integrable, P{x, v) can be normalized so that 

dv [ dxP{x,v)= [ dxP{x) = l (20) 
Jo Jo 

for all < r < 1. Note the absence of any special behavior in Eqs. (fTHj) . ()19|) at /3 = 2 or 
r = Tc. 

The probability density P(x) is shown for several values of r above and below Tc in Fig. 
2. The curves were obtained by integrating Eq. over v analytically and then performing 
the u integration numerically, using the numerical solution for P(0, u) in Fig. 1. Again there 
is no qualitative difference above and below r^- 

III. COLLISION RATE AND MOMENTS OF THE REFLECTED VELOCITY 

Unlike the distribution functions P(x, v) and P(x) considered thus far, the equilibrium 
collision rate J, defined by Eq. (jHJ, does indeed change non-analytically as r passes through 
r^. According to the asymptotic forms p4|) - (fTm) of P(0, f ), the second integral on the right 
of Eq. (jHI) converges at the upper limit for all < r < 1 and at the lower limit for (3 < 2 
but not (3 >2. Thus the boundary collision rate is finite for r > Vc and infinite for r < Tc, 
in agreement with the prediction of CSB 

The moments \v\i of the velocity just after reflection from the boundary \^ exhibit a 
closely related collapse transition. Since vP{0,v)dv is the reflected probability current in 
the velocity range v to v + dv, 

J^dvv'i+^P{0, 



1^1" = ^^^ / / , r>rc. (21) 
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The denominator in Eq. ()21|) equals the colhsion rate /, just shown to be finite for r > 
and infinite for r < r^. In the latter case, we use the regularized average 

J^dvvi+^P{0, 



kl'^ = lim ^^^TToo / , , r<rc. (22) 

From Eqs. (j2T| . (j22|l . and the asymptotic form (fT4|l . (|T5|) of P{0,v) for small f, one sees 



that all the moments \v\'^ with g > collapse at r = Tc. For r > Tc they are finite and 
nonzero, ^nd for r < rc they vanish. 

CSB [l] analyzed the case of a randomly accelerated particle, initially at x = with 
Vo > 0, moving on the half line x > with inelastic collisions at x = 0. Defining Qn{v, vo)dv 
as the probability of a velocity just after the nth refiection between v and v + dv, normalized 
so that 

POO 

/ dvQniv,vo) = l, (23) 
they calculated Qn{v,vo) and the moments 



\Vn\'^= / dvv" Qn{v,Vo) (24) 

Jo 

exactly. In the limit n oo, the gth moment diverges, independent of r, for g > |. For 
< g < |, this same quantity diverges for r > r*(g) and vanishes for r < r*(g). The critical 
parameter r*{q), given by Eq. ()15p with the replacement (3 — g + 2, decreases monotonically 
from Tc to as g increases from to |. Thus, in both the semi-infinite geometry x > 
and the finite geometry < x < 1, certain moments of the refiected velocity collapse as r 
decreases. However, since boundary collisions are less frequent in the semi-infinite geometry, 
the velocity fiuctuations are greater, and the collapse is less complete. In the semi-infinite 
case the uncoUapsed moments are infinite, the moments with g > ^ do not collapse, and for 
< g < I the critical parameter r*(g) is less than rc. 

For a particle confined to x < < 1 rather than x > 0, the recurrence relation that 
determines Qn{v,vo) is given by 



rQn+i{rv,vo) = duvG{l,v,u)Qn{u,vo) , (25) 
Jo 

Qo{v,vo) = S{v - Vo) , (26) 



as shown in the Appendix. The kernel G{l,v,u) is the same as in Eqs. ()12p and ()13p. 
Due to the property ()A2|) (see Appendix) of the kernel, the recurrence relation preserves 
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the normalization ()23p. In the hmit n — > oo, Eq. ()25|) becomes identical with the integral 
equation (fT^ for vP{0,v), suggesting that Qoo{.v,Vo) is proportional to vP{0,v). 

This proportionality could have been anticipated. In the limit n — > oo, Qn{v,VQ) is 
expected to approach the equilibrium distribution Qcqm\{v), and (5cquii('y) proportional to 
vP{0, v) follows from the interpretation of vP{0, v)dv as the reflected probability current, in 
equilibrium, in the range v to v + dv. 

The proportionality constant is fixed by the normalization ()23p . This leads to 

„ , , vP(0, v) 

Qequiil^^j = 700-1 ^TTT^ , r > r, . (27) 

For r < Tc the denominator in Eq. (|27p. which equals the collision rate /, diverges. Regu- 
larizing as in Eq. (j^ . we replace the right side of Eq. (jTTjl by limA^oQ('^^) '^)) where 

and 6'(a;) denotes the standard step function. Since dv Q{v, X) = 1, and since Q{v,X) 
vanishes in the limit A — except at f = 0+, where it diverges, 

Qeqmi(^^) = hm Q{v, A) = 6{v) , r < r, . (29) 

A — ^■O 



The distribution function QoquiK"^) collapses from (j^Tj) to as r is lowered through 



The vanishing of the moments |f I"? = 0, g > for r < Tc is consistent with the collapsed 
form 

That Qequiiiv) in Eq. (jTTjl is indeed a stationary solution of the recurrence relation 
follows directly from the integral equation (fT^ satisfied by P{0,v). That the delta function 
is a stationary solution for any r may be shown by substituting 6{u — e), e > on the 
right side of Eq. integrating over u, and then taking the limit e ^ 0. 



IV. CLOSING REMARKS 



A. Is there inelastic collapse? 

The paper of CSB ^ on inelastic collapse is almost entirely concerned with establishing 
that on the half line x > (i) the particle makes an infinite sequence of boundary collisions 
in a finite time for r < r^, and (ii) in the limit n ^ 00 the reflected velocity distribution 
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QniVjVo) and certain moments of the reflected velocity collapse as r is lowered. Our results 
for a particle in equilibrium on the finite line < a; < 1 are quite compatible with (i) and 
(ii). We question only the statement, below Eq. (19) of ^, that after undergoing an infinite 
sequence of collisions the particle remains at rest on the boundary. 

Unlike the central quantity QniVjVo) in the work of CSB, the solution P{x,v) of the 
Fokker-Planck equation ®-(IZI) provides information on both the position and velocity of 
the particle in equilibrium. The solution P{x, v) that we have obtained does not collapse 
onto the boundaries x = and x = 1 as r is lowered between 1 and 0. However, for r < r^. 



P(0, v) diverges more strongly than f ~^ in the limit f ^ 0, and this implies I = oo, |f = 0, 
g > 0, and Qeqmi{v) = S{v), via Eqs. (jS)), (EH), and (j^ . There is a collapse transition in 
the distribution of reflected velocities Qeqmi{v), but it does not involve localization of the 
particle at the boundaries. 

Why is Qeqniiiv) = 6{v) not a sufficient condition for inelastic collapse? Since the velocity 
f = on reflection from the boundary is overwhelmingly favored, doesn't the particle remain 
at the boundary? In our opinion the relevant quantity in the question of localization is not 
<5cquii(^) but the probability per unit time vP{0, v)dv for leaving the boundary with a velocity 
between v and v + dv, where 

vP{Q,v) = IQ,q,^{v) , (30) 

as in Eqs. ()27p. ()28|). If f P(0, f ) > for f > 0, the particle does not remain at the boundary. 

For r < Tc, the collision rate I is infinite, and for f > the product IQcquwiv) = 16{v) on 
the right side of Eq. is indeterminate. Whether or not vP{0,v) vanishes for f > is 
unclear from Eq. (j3Up . We have calculated P{0,v) for r < r^hj solving the Fokker-Planck 
equation. The result, as described above, is a smooth function of v with the asymptotic 
forms (fT ^ - (fT^ . The quantity vP{0,v) does not vanish for v > 0, although it does indeed 
imply Qequiiiv) = 6{v). Thus, we find that inelastic collisions do not localize the particle at 
the boundaries. 

Below we comment on two earlier results in view of these conclusions. 

B. Collision rate in simulations 

In computer simulations 0,0,0] with a discrete time step At, the boundary collision rate 
J, which can never exceed one collision per time step, is necessarily finite. In the algorithm of 
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, the root-mean-square velocity change is given by Av = {2Aty/'^. In the hmit At ^ 
the discrete dynamics approaches the continuum dynamics of Eq. and / diverges for 
r < Tc- Anton has found that the colhsion rate in his simulations scales as / ~ {At)^'^~^^^'^ , 
At ^ for r < and offered a dynamical explanation. We note that this scaling relation 
follows very simply from our results for the equilibrium distribution function P{x,v). For 
velocities |f |~Af , the simulation results are expected to deviate from the asymptotic form 
P{0,v) ~ t)~^(^') in Eqs. fll4|l . ()15|1 . Thus the boundary collision rate (jH)) in the simulations 
scales as 

/•oo 

/ ~ / dv v^-^ ~ {Avf-'' ~ (At)(2-^)/2 , At^O. (31) 

C. Persistence Exponent for r < rc 

Burkhardt and De Smedt et al. jll| have considered the probability Q(xo,fo,t) that 
a randomly accelerated particle with initial position and velocity Xq, fo, confined to the 
half-line x > and reflected inelastically at x = 0, has not yet undergone inelastic collapse 
after a time t. They predicted Q{xo, f o, t) = 1 for r > rc, and for r < Tc the power-law decay 

Q(xo,t;o;t) t^oo, (32) 

where the exponent j3 is the same as in Eqs. ()14j) , , and . In view of our conclusions 
that the particle makes an infinite number of collisions in a finite time but does not remain 
at the boundary, Q{xo,VQ,t) in Eq. (p?^ should be interpreted as the probability that after 
a time t the randomly-accelerated particle has not yet made an infinite number of boundary 



collisions. The derivations of Eq. ()32|) in 



and it is also supported by simulations 



3, HQ 



111 are compatible with this interpretation. 
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APPENDIX A: VELOCITY DISTRIBUTION ON ARRIVAL AT THE 

BOUNDARY 

The probability that a randomly accelerated particle with initial position a; = and 
initial velocity u > 0, moving on the half line x > 0, arrives with speed between v and v + dv 
on its first return to a; = is given by vGq{v, u)dv, where 

Goiv,u) = — Al 

This result, due to McKean P|, was also obtained independently by CSB Q]. 

The quantity G{l,v,u) in Eq. (fT^. derived by BFG [4j, extends this result to the the 
finite interval < x < 1. The probability that a randomly accelerated particle which leaves 
a; = with velocity m > has speed between v and v + dv the next time it reaches either 
boundary is given by vG{l,v,u)dv, where the first and second terms on the right side of 
Eq. (fTSjl correspond to arrival at x = and x = 1, respectively. Like Gq{v,u) in Eq. ()A1|1 . 
G{l,v,u) satisfies the normalization condition 

dvvG{l,v,u) = 1. (A2) 







Integral equation (jl2|) for P{0,v) follows directly from the interpretation of G{l,v,u) 
in the preceding paragraph and the stationarity, in equilibrium, of the reflected current 
vP{0,v)dv between v and v + dv. Another consequence is the recurrence relation (j25|) for 
the probability distribution Qn{v,VQ) of the speed with which the particle rebounds after 
the nth boundary collision. Solving Eqs. ()25|1 . (f^H|) with Gq{v,u) in Eq. (jAlll in place of 



G{l,v,u), CSB [li calculated Qniv,Vo) exactly for motion on the half hue x > 0. 
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FIG. 1: Double-logarithmic plot (base 10) of P{0,v) for r = 0.01 (dotted curve), r = 0.1 (solid 
curve), r = 0.2 (dashed curve), and r = 0.5 (dot-dashed curve). The curves are normalized 
according to Eq. ((^OJ. 
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FIG. 2: Double-logarithmic plot (base 10) of P{x) for several values of r. The curves are normalized 
according to Eq. (|2()|) . 



